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We consider a one-dimensional spin chain for which the ground state is the cluster state, capable 
of functioning as a quantum computational wire when subjected to local adaptive measurements of 
individual qubits, and investigate the robustness of this property to local and coupled (Ising-type) 
perturbations. We investigate the ground state both by identifying suitable correlation functions as 
order parameters, as well as numerically using a variational method based on matrix product states. 
We find that the model retains an infinite localizable entanglement length for Ising and local fields 
up to a quantum phase transition, but that the resulting entangled state is not simply characterized 
by a Pauli correction based on the measurement results. 
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I. INTRODUCTION 

Measurement-based quantum computation 
(MBQC) [lj has recently emerged as an alternative 
model for quantum computation to the standard circuit 
model [2 . In MBQC, computation proceeds via local 
(single-qubit) adaptive measurements on a fixed fiducial 
state of a quantum many-body system, for which the 
cluster state [TJ [3] is the canonical example. 

The MBQC model of quantum computation is partic- 
ularly useful for investigating the physical requirements 
for quantum computing, as it becomes possible to pose 
questions about universality and fault-tolerance in terms 
of properties of quantum states [1 IS] l6[ l7[ IS] l9l flT) l IPT ] 
1121 113] . For example, the universality of a given quantum 
state for MBQC may be determined by assessing the fi- 
delity and range of a universal quantum gate set, which 
in turn is quantified by viewing gates as resource states 
for gate teleportation [3J |T3] prepared via local measure- 
ments [15] . 

However, despite this powerful framework, there are 
relatively few known examples of resources states (or dis- 
tinct classes of states) that allow for MBQC [HJ [71 [S]. 
A promising avenue for identifying properties of states 
that allow for MBQC is by investigating ground or low- 
temperature thermal states of a coupled quantum many- 
body system 0[ini[T7]. In particular, one can construct 
model Hamiltonians for which the ground state is univer- 
sal for MBQC (say, the cluster state on some appropriate 
lattice) and then investigate the robustness of this prop- 
erty to local perturbations. Progress has been made in 
this direction for thermal or local perturbations of the 
cluster state PlfTUlfTTj. 

In this paper, we investigate a ID chain of qubits 
for which the cluster state is the ground state, and in- 
vestigate the robustness of its computational power to 
both local and coupled (Ising-type) perturbations. (For a 
ID chain, "computational power" is restricted to single- 
qubit unitary evolution.) Consistent with previous re- 



sults [10] lllj. we identify a robust phase for which every 
ground state can allow for quantum information to be 
transferred using local measurements. We investigate the 
usefulness for states in this phase to serve as "quantum 
computational wires" |18j . which function as primitives 
for MBQC. 

A. Quantum computational wires 

Within MBQC, the simplest primitive is the ability to 
move information (i.e., teleport) along one-dimensional 
channels, with single-qubit unitaries determined by the 
choice of measurements. States with this property are 
known as quantum computational wires [18 . A wire in 
MBQC consists of two parts: (i) creating a maximally- 
entangled state between two distant points via local mea- 
surements; and (ii) identifying the correct "bi-product" 
unitary based on the measurement result (either to im- 
plement the identity gate, or some more general single- 
qubit unitary gate). The first property is characterized 
by localizable entanglement [19] , which is the maximum 
average entanglement that can be localized on two sites 
through local measurements on all other sites. For sys- 
tems where the localizable entanglement falls of exponen- 
tially, the entanglement length £e is defined through 

E L ~ e~ n/iE when n > 1 , (1) 

with n the separation between the sites; states with fi- 
nite localizable EL are not directly useable as a quantum 
computational wire. In contrast, systems with infinite lo- 
calizable EL allow for teleportation over arbitrary length 
scales. One of the characteristics of the cluster state is 
a diverging localizable EL. (We note, however, that di- 
verging EL is in itself not necessary nor sufficient for the 
state to be useful for MBQC 0.) 

In general, identifying an optimal measurement basis 
for localizing entanglement can be a challenge. For the 
cluster state on any lattice, such an optimal measure- 
ment strategy is known: measure Z (i.e., measure in the 
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eigenbasis of a 2 ) on all qubits except those on a line con- 
necting the two desired qubits, thus effectively making a 
ID cluster state, and then measure X on all intermediate 
qubits along the remaining line. This measurement se- 
quence will concentrate a maximally-entangled state on 
the two remaining qubits. 

To transform the resulting maximally-entangled state 
into a particular one (say, the two-qubit cluster state), a 
correction unitary on one of the qubits based on the mea- 
surement results is then applied. If these corrections are 
Pauli operations, or in general if they close to form a fi- 
nite subgroup, then the bi-product operators do not need 
to be performed but instead one can use a classical com- 
puter to keep track of them. In general, however, even 
if the localizable EL is infinite, the bi-product operators 
may not close in this way. That is, there is a distinction 
between infinite localizable EL and the ability to function 
as a quantum computational wire. 

This paper is structured as follows. In Sec. [TTJ we 
introduce our Hamiltonian with local and coupled per- 
turbations, and investigate its phase diagram. Here, we 
identify a cluster phase connected to the cluster state 
(with zero perturbations). We consider using correlation 
functions as for the cluster state to quantify the identity 
gate within the cluster phase in Sec. |III| In Sec. |IV] we 
use numerical methods to assess the localizable entan- 
glement within our model, and in Sec. [V] we investigate 
the usefulness of ground states in the cluster phase to 
serve as a quantum computational wire 
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In Sec 

turn our attention to the use of local measurements to 
disentangle the two halves of the chain, in analogy to the 
Z measurement in cluster state MBQC, and finish with 
conclusions in Sec. |VIl[ The appendix provides details on 
the Jordan- Wigner transformation. 



II. CLUSTER HAMILTONIAN WITH 
ANOMALOUS TERMS 

A. The cluster state and cluster Hamiltonian 

Consider a graph (such as a lattice) containing N ver- 
tices, with a qubit placed on each vertex. The cluster 
state on this graph can be constructed by first prepar- 
ing each qubit in the state |+), the +1 eigenstate of the 
a x Pauli spin operator, and then applying a controlled 
sign operator U = exp(wr|l)(l| ® on every pair of 

qubits on vertices connected by an edge. The resulting 
cluster state is compactly described in terms of stabiliz- 
ers. That is, at every vertex /i = 1, . . . , N one associates 
an operator 



K„ 



n 



(2) 



where v ~ \i indicates all vertices v that are connected 
to fi. In a graph of degree %, this operator acts on \ + 1 
vertices. The cluster state is the unique state \C) that 



*gc> = |c>, v M . 



(3) 



An alternate approach to prepare a cluster state would 
be to view the qubits on this graph as an interacting 
quantum many-body system, and to find some Hamilto- 
nian He for which the cluster state is the unique ground 
state. Provided the model has a sufficiently large gap, 
preparation of the state would amount to cooling the 
system down to (or near) the ground state. One such 
Hamiltonian that has the cluster state as its ground state 
is (minus) the sum of the stabilizers at every site \i [4] , 



z2 K » 



(4) 



We have chosen units such that the energy scale is fixed, 
and the gap in the model is 2 between the unique ground 
state and a TV-fold degenerate first excited state. Al- 
though the terms is this Hamiltonian are gener- 
ally many-body, such interactions can occur as the low- 
energy behavior of a more natural two-body Hamilto- 
nian [Tni [T7]. We note that the cluster state on a line - 
the ground state of the Hamiltonian Q - has an infinite 
localizable EL and can serve as a quantum computational 
wire. This Hamiltonian on a line can be realized in an 
optical lattice despite its three-body interactions [20 . 



B. Cluster Hamiltonian with anomalous terms 

Consider the cluster state Hamiltonian Q in one di- 
mension with two types of additional terms: local fields 
and couplings, 



H = Hr 



(5) 



with J > 0. Unless otherwise specified, we will con- 
sider periodic boundary conditions, such that (Jn+i = 
Throughout the article we will restrict our attention to 

By = 0. 

(Note that an Ising interaction as in the last term of 
Eq. ([5]) can be used in a time modulated fashion to con- 
struct a cluster state [3]. However, we will study this 
term only as a constant perturbation of the Hamilto- 
nian.) 

A number of results are known for specific cases in- 
volving a single local term and J = 0: 

Local B z field: In one dimension, the system is fun- 
damentally unstable with the addition of a local z field, 
i.e., the EL becomes finite for any non-zero B z . In two- 
or higher-dimensional lattices, however, it can be shown 
using techniques from percolation theory that the model 
exhibits a transition from a finite region of parameter 
space B z < B z rlt wherein the ground state has infinite 
localizable EL. Because the localizable EL is zero in the 
limit B z — > 00, there is a transition in the localizable EL 
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for two- or higher-dimensional lattices even though the 
underlying model does not exhibit any quantum phase 
transition [TTj . 

Local B x field: In one dimension, the system exhibits 
a single quantum phase transition at B x — 1 separat- 
ing the "cluster phase" and a separable phase. Ground 
states in the cluster phase are characterized by infinite 
localizable EL [TCJl I2TJ] , and the system serves as a quan- 
tum computational wire at all length scales (with pre- 
cisely the same measurement sequence and corrections 
as for the cluster state) albeit with lower fidelity for all 
B x < 1. The performance as a quantum computational 
wire is therefore a robust property of this system in the 
presence of a perturbing B x field. 

With these prior results, we now analyze the full phase 
space of the Hamiltonian Consider the action on 

this model of the unitary transformation U = FJ ^ W 
that applies the controlled phase gate E/ M „ on all pairs of 
adjacent qubits. This unitary maps He to 

Uif c Ut = -5>*, (6) 

for which the ground state is separable, U|C) = |+)' 8Ar . 
In general, the transformation leaves a z operators invari- 
ant, while <t x maps to 

U^Ut =a *Y[ot = K M . (7) 

This transformation allows us to determine some prop- 
erties of the model First, consider the model with 
B = 0. Under this unitary mapping, the model is dual to 
the ordinary transverse-field Ising model, which is com- 
pletely solved in ID [21], and has a single quantum phase 
transition at J — 1. For our model ([5]) in the large J limit 
with B = 0, the ground state approaches a GHZ state 

|G#Z) = ~(|00---0) + |ll---l)) . (8) 

We denote the phase for J > 1 the Ising phase. 

From the local B z results, we do not expect the proper- 
ties of the cluster state to survive for any B z > 0, and so 
we now consider the restricted model with B z = 0. This 
model can be subjected to a Jordan- Wigner transform as 
shown in Appendix [A] which easily allows identification 
of critical points in the model, and hence diverging corre- 
lation and entanglement lengths. This restricted model 
has three clearly identifiable phases: the cluster phase 
as B X ,J — > 0, an Ising phase as J > max.(l,B x ) and a 
separable phase as B x 3> max(l, J). The model is uni- 
tarily dual under the transformation U. Thus, for any 
critical point at (J,B X ), there is another critical point at 
(-jj-, -g-). Thus, considering the critical line Cc connect- 
ing the known critical points at (0, 1) and (1, 0), this dual- 
ity reveals another critical line Cq from (0,1) to (oo, oo). 
Parameterizing the lines with a parameter r S [0, 1], we 




FIG. 1: A sketch of the phase diagram for the Hamiltonian in 
Eq. |5| with B z — 0. The lines Co and Cc constitute quantum 
phase transitions between the three indicated phases. (All 
quantities shown in figures in this paper, including this one, 
are dimensionless.) 



have 

C-fc/M) fe-^,^), (9) 

for some unknown function /(r), for which /(0) = 1 and 
/(l) = 0. In the limit B X ,J 1, the cluster term in 
the Hamiltonian becomes unimportant, and the model 
is a simple Ising model for which the critical point must 
be at B x = J. Under the above parameterization we 
must have, assuming /(r) is continuous, and invoking 
l'Hopital's rule, 

r 1 t - 1 1 

lim — — - — — — lim = lim = , (10) 

f( T ) f( T ) T ^ f( T ) T ^ X d f/ dT 

which implies ^ = oo. Hence, the critical line Cc must 
be convex (at least close to J = 1). The numerical results 
indicate that this convergence is very slow, and to a close 
approximation we have f(r) = 1 — r, and hence Cq = 
(t', 1+t'), t' £ [0, oo]. The three phases are separated by 
Ising type critical lines corresponding to a conformal field 
theory with central charge c= 1/2, and a corresponding 
entanglement signature 22J. 

In the phase diagram of Fig. [T] for B z = 0, at the origin 
(J = B x = 0) the ground state is the cluster state. From 
the results of (TUJ |20j, we know that ground states on 
the line J = 0, < B x < 1 have infinite localizable EL, 
and in fact allow for long-ranged single-qubit gates using 
the same measurement sequence and Pauli corrections as 
for the cluster state. We now turn our attention to the 
question: is any ground state within the cluster phase of 
Fig.[l]also useful as such a quantum computational wire? 
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III. CORRELATION FUNCTIONS FOR 
LOCALIZABLE ENTANGLEMENT AND 
QUANTUM COMPUTATIONAL WIRES 

In special cases, the localizable entanglement of a state 
can be characterized (or, at least, lower bounded) by a 
correlation function. For example, with certain measure- 
ment sequences, one can obtain post-measurement quan- 
tities from pre-measurement expectation values to quan- 
tify the localizable entanglement [H)l |2"51 |2"4"] or the per- 
formance of quantum gates in MBQC [lOj [15]. That is, 
the expectation values of string-like operators can serve 
as order parameters to identify the hidden correlations 
corresponding to localizable entanglement and the abil- 
ity to function as a quantum computational wire. 



A. Localizable entanglement for B z = 

For B z = 0, the symmetries of our model allow us 
to determine the optimal measurement sequence and 
to derive an expression for the localizable entangle- 
ment [19J HA [24]. Our model ^ with B z = is in- 
variant under rotations of all spins by tt about the x 
axis, transforming <r z — > —<r z . It is not invariant under 
similar rotations about the y and z axes; therefore, the 
optimal measurements for localizing entanglement on a 
finite chain between the end qubits 1 and N arc measure- 
ments in the X-basis QH H3J HQ. With this basis > thc 
localizable entanglement for this finite chain is given by 
the string correlation function 

^ = K(nf= 2 V)o. (ii) 

(We note that, for a suitable choice of boundary condi- 
tions on this finite chain, we could ensure that the ground 
state is a +1 eigenstate of P x — JI^Li a j an d simplify 
this correlation function further. However, as we will 
be interested in open chains or periodic boundary condi- 
tions, we will not do so.) 

This correlation function quantifies thc localizable en- 
tanglement, but does not completely characterize the 
form of the resulting maximally entangled state. We now 
turn to a complete set of correlation functions for this, 
which quantify the state's use as a quantum computa- 
tional wire. 



B. Correlation functions for quantum 
computational wires 

In this section, we briefly summarize the main result of 
Ref. [15]. Consider a lattice of qubits prepared in some 
initial pure state \ipo}- Singling out two qubits, a and 
b, we consider a measurement sequence on the remain- 
ing qubits in the lattice that localizes entanglement on 
a and b. Let m label the measurement outcomes, and 



P m be the corresponding projector. Following the mea- 
surements, a correction unitary U m conditional on m is 
applied to qubit b. Averaged over all possible measure- 
ment outcomes, the resulting two-qubit state is 

Pat = Y. m u m PM(MP m ul . (12) 

Equivalently, we can characterize this final state using 
the set of expectation values of bipartite Pauli operators 
a a® °6 j h 3 — It x tVt z on qubits o and 6, as 

= E m (HPm<(4) m P m m , (13) 

where {<J 3 b ) m = U^alUm, and where (-)-p denotes the ex- 
pectation value in thc final post-measurement two-qubit 
state. The set of such correlation functions for all pairs 
of Pauli operators will completely specify the two-qubit 
state. 

We now restrict our attention to measurement se- 
quences for which there exists a string of operators S 
acting on some set of the measured qubits which is inde- 
pendent of the measurement outcomes, and an operator 
tI on b that is also independent of m, such that 

<K) m Pm = P m <Srl . (14) 

For example, in the cluster state model of MBQC [3J, 
a universal gate set is known that satisfies this prop- 
erty |15j . For such measurement sequences, using the 
projector properties P^ = P m and J2 m P ™ = 1 g ives 

(aiai) v = (Vdk^'lVo) ■ (15) 

Thus we can relate the two-qubit state prepared after the 
sequence of measurements to a correlation function of the 
original state |^o) prior to measurements. That is, the 
correlation functions characterize the jwsi-measurement 
two-qubit state using expectation values of strings of op- 
erators on the pre-measurement state. 

It is critical to this development that one can identify 
such a string of operators S. The measurement sequence 
for localizing entanglement in thc cluster state provides 
the canonical example. For a iV-qubit chain with N even, 
using the same measurement sequence as if localizing en- 
tanglement in the cluster state |C) (measure Z on qubits 
1 and N, and X on all qubits 3 • • • N — 2) with the same 
Pauli corrections as for the cluster state, the averaged 
state on qubits a = 2 and b = N — 1 afterwards has the 
following non-zero correlation functions: 

(aX)p = (^kf (nfii~H+iHl<?> (16a) 

(o* a oi) v = (ci*f(nfi?~H-H-ii c '> ( 16b ) 

{alal) v = (Ck^nf^^^lC) (16c) 

Note that a two qubit cluster state has expectation val- 
ues {o x a a z ) — (cr z a x } — {a v a a\) = 1, which means that 
the above expectation values are maximal for the cluster 
state. 
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FIG. 2: Expectation values (a^a^)-p and (a^a^)-p for B x = 
B z = 0. The former is shown for N = 100, 200 while {a y a <J v b )v 
is shown for TV = 200 only, as this is not size dependent, 
and indistinguishable from TV = 100. The expectation value 
(pl&b)v is identical to (o^oDv- 



C. Correlation functions for B z = 

For our model (|5| with B z = 0, the correlation func- 
tions of Eq. ( fl6| can be calculated analytically using the 
exact solution obtained via the Jordan- Wigner transfor- 
mation. In the cluster phase, of the three correlation 



functions given in Eq. (16 1, we find that only the expec- 
tation value (u v a a\)-p of Eq. (16c) remains long-ranged in 



the thermodynamic limit. This expectation value takes 
the general form of a transverse-field Ising order parame- 
ter. The other expectation values, {o z a o^)v and (a^a^)-p, 
remain nonzero only for J = 0. (As shown in JlOj . for 
J = they remain long-ranged for B x < 1.) Fig.[2]shows 
the behaviour of these expectation values for the ground 
state of the model with B = 0. 

The long-ranged behaviour of (a v a a\)-p in the cluster 
phase serves as a useful order parameter, and is related 
to the expression for the finite-chain localizablc entan- 
glement of Eq. (111. However, the fact that the other 



correlation functions of Eq. ( 16 1 are zero for J > shows 



that the standard cluster-state Pauli corrections based 
on the measurement results are not suitable, i.e., the re- 
sulting bipartite state is not the two-qubit cluster state. 
We now turn to numerical methods to investigate the ex- 
act form of the resulting bipartite entangled state in the 
cluster phase. 



IV. NUMERICAL CALCULATIONS USING 
MATRIX PRODUCT STATES 

Our analytic approach to investigating the cluster 
phase has identified an appropriate order parameter for 
the localizable entanglement, realized as the expectation 
value of a string operator, but has not revealed the form 



of the resulting bipartite entangled state necessary for use 
as a quantum computational wire. We now use numerical 
methods based on Matrix Product States to investigate 
the localizable EL in this phase. 



A. Matrix Product States 

Matrix Product States (MPS) have emerged as the nat- 
ural language to describe ID systems with limited long- 
ranged entanglement, i.e., off criticality |25j . (There are 
examples of MPS critical points where the states sub- 
scribe to an exact MPS representation [26], but general 
quantum critical points will be associated with a diverg- 
ing entanglement.) The well-established density matrix 
renormalization group (DMRG) scheme [37] has subse- 
quently been shown to be an iterative minimization over 
MPS. Such states are also well suited to describe lo- 
cal measurements, and hence naturally suited to cluster 
states and measurement-based quantum computation in 
general [6] [7J . Even when a full analytic solution is avail- 
able through the Jordan- Wigner transformation, we still 
find that an MPS representation is more useful for this 
reason. 

Given a sufficiently-large bond dimension D, any state 
can be written as a MPS of the form 

d 

|V.) = J2 T±(A?-.A' ] ?)\8 1 -..8„), (17) 

8f-S N =l 

where d is the local Hilbert space dimension (e.g., d = 2 
for qubits). Each site \x is associated with d different 
D x D matrices ^4^, s = 1, 
to either of 



, d, which are normalized 



5>*(^)t=l, (18) 



for all (i. For a translationally invariant state, we have 
A"^ = A s independent of fj,. 

We now consider MPS descriptions of the ground state 
of our Hamiltonian Q in one dimension. The cluster 
state is described by the simplest nontrivial case, with 
D = 2 and 



.4 U 







A 1 



(19) 



With J = B x — but B z > 0, one still has a D = 2 rep- 
resentation. Any MPS representation can be obtained if 
one knows a sequential procedure to produce the state 
such that two neighbouring states are a result of opera- 
tion by a two-qubit operator U when all qubits initially 
are in the state |0) [25]. In the cluster state U — U(H®1) 
where H is the Hadamard transform and U is the con- 
trolled phase gate. For the B z > case, the controlled 
phase gate commutes with the extra term in the Hamil- 
tonian, so the resulting ground state can be written as 



6 



\C Z ) = U\6)® N where |0) = cos6>|+) + sin6»|— ) and 



tan ( 



(20) 



Hence, to obtain the MPS one can use the same sequen- 
tial generation scheme, replacing the Hadamard by 



H 9 = 



V+ V- 

n- -v+ 



with r)± = (cos(9) ± siu(0))/v% which results in a MPS 
representation 



AO _ V+ V- 
• 1 



A 1 





-v- 



(21) 



All other points in the phase diagram have only approx- 
imate solutions for any D < 2 N ' 2 . 

Given a measurement sequence over the qubits 1Z with 
outcomes Si,i € 1Z, the resulting state after measurement 
with the measured qubits traced out is 



Tr (Al 1 



(22) 



where the sum is over all i 1Z. A general measurement 
on qubits is given by a direction on the Bloch sphere. As 
we have restricted our attention to B y = 0, the ground 
state of our Hamiltonian will always have real coefficients, 
and we can restrict our measurements to directions in the 
x — z plane, given by a direction £. We define £ = (tt/4) 
to correspond to an Z (X) measurement. 

An optimal measurement basis which localizes the 
maximum entanglement is given by the angle that max- 
imizes J2i I det A j\ [29] . Applying this result to states of 



the form of Eq. (21 1, the MPS representation in a tilted 
basis is 



A { ^ = cos £A + sin £Ai , Af ' = - sin £_A + cos £Ai 



(?) 



and we get 



(23) 



Ei det ^ (c) i 



|2sin(2£)r7 + 77_| = | sin(2£) cos (20) | 



Thus, an optimal measurement basis for localizing en- 
tanglement in this state is for £ = 7r/4, i.e., the X basis, 
independent of the magnetic field B z . 



B. Variational method 

The problem of finding a MPS representation for the 
ground state of a given Hamiltonian is an NP-completc 
problem [3D]. However, variation methods such as 
DMRG work well in practice for a wide variety of Hamil- 
tonians. In addition, the MPS representation makes 
it easy to compute expectation values, entanglement, 
and other physical quantities, and make comparisons to 



known values where such exist. The numerical successive 
minimization is described in other works [19 [ I28 ] I5T]. but 
a short description of what we denote the Variational 
Matrix Product State (VMPS) procedure follows. Start- 
ing with a random MPS representation and correspond- 
ing energy, one fixes all matrices except one, and mini- 
mizes the energy with respect to the single matrix. This 
minimization amounts to a simple generalized eigenvalue 
problem. After normalization of the matrix, one moves 
to the neighboring site, and then repeatedly sweeps the 
lattice back and forth until convergence is reached. We 
stop the iteration after a sufficient number of full sweeps 
plus half way back in the lattice where we pick out the 
two matrices needed, and use these as a representation for 
a large chain with periodic boundary conditions. Thus 
we avoid boundary effects, and we get a convergent rep- 
resentation. 

The performance of the method is naturally dependent 
on the initial choice of representation, and it is useful to 
repeat the procedure with different choices to find a con- 
sistent solution. In general, energy and local expectation 
values converge quickly and independently of the initial 
state to the same value. However, even very subtle differ- 
ences in these quantities can amount to a large differences 
in terms of entanglement quantities. 

We can assess the accuracy of our VMPS for B z = 0, 
by comparison with the analytic solution available via the 
Jordan- Wigner transform. With the exact solution, we 
can easily compute any n-partite reduced entropy S n = 
S(p n ), where p n — Tr n+ i...N\ip)(ip\ is the reduced state 
of n neighboring qubits for a system in the pure state 
and S(p) is the von Neumann entropy. The singlc- 
qubit reduced entropy Si will not distinguish the cluster 
and Ising phases, as both will have Si = 1, so we use 
the bipartite reduced entropy as an indicator. This 
entropy separates the three conventional phases of S2 = 2 
(cluster), 52 = 1 (Ising) and S2 — (separable). A state 
that is very close to the true ground state in terms of the 
energy and expectation values can still be very different 
in terms of 52, so it is necessary to run the iteration 
several times with different initial condition to obtain a 
state that reflects the ground state in this respect. 

We consider the line B z = 0, J = 0.5, where we ex- 
pect two Ising class phase transitions at B x w 0.5, 1.5. 
The bipartite reduced entropy for a number of runs with 
random initial conditions are shown in Fig. [3] where it 
is clear how the entanglement properties fail for ground 
states in the Ising phase (which are GHZ-like), while be- 
ing relatively good for the other two phases. The local 
expectation value (a*) is however well represented for all 
phases. These results also show that 52 is upper bounded 
by the true value, and the best representation can be 
chosen as those with the highest 52 even when the exact 
value is not known. 
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FIG. 3: The bipartite reduced entropy S2 for J = 0.5, B z = 
computed exactly through the Jordan- Wigner transform (full 
line) and through 40 runs of VMPS with 6 sweeps. Note 
that the cluster and separable phases compute the entangle- 
ment properties very precisely, especially away from the phase 
transitions, while in the Ising phase the results are much less 
reliable. Also, the algorithm inevitably underestimates the 
true value of 52 • The inset shows the local expectation value 
(ct^), where error bars are too small to be evident. (N=200, 
D=8) 



FIG. 4: The localizable entanglement for a N = 60 model 
with 5 = 0. MPSs used are D = 4 (open circles) and D — 8 
(solid circles). The gray line is a guide to the eye, follow- 
ing the function (1 — J 2 ) ' 36 , which is consistent with the 
expected behaviour of an second order phase transition or- 
der parameter. The error bars indicate the sensitivity of the 
initial conditions in the VMPS scheme. The inset shows a 
similar plot for B x = .5 and D = 8. Note that the perfor- 
mance of the VMPS is noticeably worse as J approaches the 
phase transition, though the data are still conclusive. 



C. Numerical estimation of localizable 
entanglement 

We now seek find an approximate MPS representation 
for the ground state of our model and with it, to com- 
pute the localizable entanglement. The latter step can 
be accomplished by a Monte-Carlo scheme [TS] to sam- 
ple over entanglement by a weighted random walk in the 
probability space. Given that two substantial numerical 
steps are needed in this method to obtain the localiz- 
able entanglement, and in particular as a tendency of the 
MPS construction to occasionally yield incorrect repre- 
sentations of the ground state, the values obtained has 
an intrinsic error that we have indicated wherever ap- 
plicable. Specifically, in the Ising phase, the method is 
unstable and picks out a ground state close to a product 
state rather than the GHZ state. However, we are in- 
terested in what happens in the cluster phase, and that 
problem is therefore not substantial. Importantly, the 
technique seems to be particularly well behaved in the 
cluster phase. 

We use periodic boundary conditions (PBC) through- 
out, as this makes the computation much easier and, due 
to the resulting translation-invariance, reduces storage to 
only two matrices as opposed to 27V matrices in the open 
boundary condition case. As we show in Sec. |VI[ the Z 
measurements do not perfectly disentangle the chain for 
J > 0, and one might worry that the results are depen- 
dent on the boundary conditions as there are two direc- 
tions in the lattice that the entanglement might propa- 
gate. However, our data show that this does not happen. 



Specifically, one can do n z subsequent Z measurements, 
and for n z > 1 , the resulting values for the entanglement 
is independent of n z for all J. Hence we are justified to 
use PBC, and the error in doing so is much smaller than 
the statistical noise of our methods, and we use n z = 2 
in the following. 

To estimate the localizable entanglement we must 
choose a specific measurement protocol, and the results 



of Sec. Ill reveal that the optimal basis for measuring 
the intermediate qubits is the X basis. We measure the 
qubits beyond the endpoints in the Z basis as for the 
cluster state. 

Our numerics confirm our expectation that the localiz- 
able EL becomes finite for any B z > (as can be shown 
analytically at J = B x = using the MPS representation 
of Eq. 21 1. For the remainder of this paper, we restrict 
to B z — 0. We have computed the localizable entangle- 
ment along two lines in the phase diagram of Fig. [T| (i) 
B x = B z = and (ii) B x = .5 and B z = 0. The (i) line 
will have an Ising class quantum phase transitions sepa- 
rating the cluster and Ising phases, while the (ii) line will 
have a phase transition at roughly J — 0.5 between the 
same phases. 

For the model with B = 0, the data clearly shows 
an infinite EL for J < 1, and a significant localizable 
entanglement, as shown in Fig. [4] Similar results are 
found along the line (ii), with B x = 0.5. 

The sensitivity on the initial conditions in the VMPS 
algorithm is far larger than the numerical errors of the 
Monte Carlo procedure. However, the latter is well suited 
to detect polynomial versus exponential decay, and per- 
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forms consistently in this respect for different initial con- 
ditions even though the actual entanglement values may 
have substantial variation. Hence, we are able to sepa- 
rate an infinite EL phase (the cluster phase) and a finite 
EL phase. 



V. CHARACTERIZING ENTANGLEMENT IN 
THE CLUSTER PHASE 

We have shown both analytically and numerically that 
the localizable EL is diverging in the cluster phase, and 
the numerics also provide a mechanism to analyze the 
form of the resulting two-qubit entangled state. For 
J > 0, the resulting two-qubit state is not identical to 
the form of the two-qubit cluster state even after the 
Pauli corrections. In this section, we explore the form of 
the post-measurement state and its dependence on the 
measurement results. We focus specifically on B = 0, 
but our results extend directly to B x > 0. 

As we are performing the identical measurement se- 
quence for localizing entanglement in the cluster state, 
it is natural to apply the same Pauli corrections to the 
final two-qubit state depending on the parity of the even 
and odd X-measurements, as well as the boundary Z- 
mcasurcmcnts. The correlation functions ( 16 1 which in- 



corporate these corrections show that the resulting state 
does not take a fixed form. With the above Pauli cor- 
rections, the sampled states in the Monte Carlo anal- 
ysis with near-maximal entanglement are characterized 
by {a v a u\)-p = 1, as expected from our analysis using the 
Jordan- Wigner transformation, as well as the relations 



K°l)v, (25) 



and 



z _x\ 2 
J a (J blV 



J J\2 
7 a a b)V 



= 1 . 



(26) 



(Sampled states that are not near-maximally entangled 
do not fall into this class.) The conditions (251 imply 



that the maximally entangled state are described by 

|$) = cos$|C , oo)+sin$|C 11 ) ) (27) 
where we have defined 



\ c u) = (°f )' K) J |c) i,i = o,i. 



(28) 



This form of the states is consistent with the long- 
ranged behaviour of (o^er^)-p, but reveals an additional 
y-rotation by an angle 

Investigating the distribution of angles $ for states 
sampled from the Monte Carlo analysis, we find that 
these angles do not take fixed multiples of ir/2. This state 
may be corrected into the cluster state by a y-rotation of 
qubit b by an angle — $, but not by a Pauli correction. 
The angle $ is a function of the measurement results, 
but we have been unable to determine this dependence. 




FIG. 5: The amplitude A(J) for the sinusoidal behavior of 
(xz) after a projection P and SO(2) correction on one qubit 
is done on the state. Closed circles are for N = 50 while open 
circles are N = 102. The inset shows the probability p yy of 
success for the measurement. Note that in the latter case the 
two system sizes are practically indistinguishable. 



We model <I> as sampled from a probability distribution 
Vj{<&), satisfying X> ($) = <$($). We now show that this 
distribution does not possess any bias away from zero; if 
that were the case, one could improve the fidelity of the 
post-measurement state with the two-qubit cluster state 
by performing a correcting rotation. For a given J, the 
expectation value (a^a^)-p will take the form 



biv 



A(J)cos(e(J)), 



(29) 



where ,4(0) = 1, f (0) = 0. The phase £(J) thereby deter- 
mines the bias, and the amplitude A gives the magnitude 
of the expectation value. 

We can now investigate the behaviour of our localized 
entangled state as we approach the phase transition at 
J = 1. We restrict our sampled states to maximally- 
entangled ones by projecting onto the subspace of states 
|$), i.e. using a projector P = |Coo)(C o| + |di) (di |. 
The probability of this projection viewed as a measure- 
ment is denoted p yy ; we note that this probability is large 
for all J up to the phase transition as seen in the inset 
of Fig. 

On states following the projection, we investigate the 
average angle. We find that the phase £( J) is an indicator 
of a size-dependent transition, changing rapidly from zero 
to it/2, with a behavior closely approximated by 

f(J) - \ (t & nh{K(N)[J ~ V (N)}} + 1) , (30) 

where the coefficients K and r\ are size dependent. The 
parameter r\ is a good marker for the transition, and as 
shown in Fig. [H] it follows a power law ?y ~ iV -1 / 3 , while 
A increases non-universally with N. Up to this transi- 
tion, there is no bias, and the state approximates a cluster 
state. Beyond this transition, the distribution T>j(Q) is 
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FIG. 6: The phase £(J) is shown for AT = 50 (closed circles) 
and N = 100 (open circles). The error bars refer to the fitting 
to the sinusoidal curve. Full lines are best fits to tanh, while 
the inset shows the scaling law for the transition parameter 
r\ ~ N~ 1 ^ 3 on double log scale. 



almost uniform, meaning that there is no rotation inde- 
pendent of the measurement results that will increase the 
fidelity with the two-qubit cluster state. 

In summary, while the localizablc EL remains infinite 
throughout the cluster phase, the resulting two-qubit en- 
tangled state cannot be deterministically transformed to 
the fiducial two-qubit cluster state using Pauli correc- 
tions for any J > 0. Further work is required to deter- 
mine the non-Pauli correction angle (about the y-axis) 
based on the measurement results. Potentially, this non- 
Pauli rotation could be useful for developing quantum 
computational wires |18j that perform non-trivial singlc- 
qubit gates. However, without characterizing this rota- 
tion in terms of the measurement results, such states can- 
not be used directly as quantum computational wires. 



VI. DISENTANGLING MEASUREMENTS 

Finally, in addition to localizable entanglement, we 
consider another property of quantum states that is use- 
ful (though possibly not necessary) for measurement- 
based quantum computation. The cluster state possesses 
the useful property that a Z measurement on any qubit 
"removes" that qubit and leaves the remaining qubits in 
a cluster state (up to a Pauli correction dependent on 
the measurement result). On a 1-D chain, a Z measure- 
ment on a qubit will then disentangle the two halfs of the 
remaining chain. This property is shared by our model 
Hamiltonian for any B z as long as B x = J = [1 1 . If we 
turn instead to the model with B = but J > 0, we find 
that on an N — 3 open chain, any measurement on the 
middle qubit will result in some remaining entanglement 
on the two end qubits, and therefore be an inadequate 
disentangling measurement. 




FIG. 7: For a chain of 4 qubits, the measurement angle that 
disentangles the two extreme qubits is shown as a function 
of J in the main panel. The left inset shows the correspond- 
ing outcome probabilities for the three possible measurement 
outcomes. The right inset shows the remaining entanglement 
multiplied with the probability of the outcome between the 
end qubits if no correction is done. In this case the outcomes 
s = ±1 are identical, so there are only two distinct cases, and 
the average entanglement is shown as a grey line. 

In this section, we consider using local measurements 
on pairs of qubits that optimally disentangle the remain- 
ing halves of a 1-D chain with B — 0. (To be clear, we 
still consider single-qubit measurements, applied to two 
neighbouring qubits.) First, we consider both qubits to 
be measured at an angle i9 in the x—z plane. The out- 
comes s are labelled according to the total spin in the de- 
sired direction which is either of s = {— 1,0,1}. There are 
two outcomes with total spin zero, but these lead to iden- 
tical states due to the reflection symmetry of the model 
and are therefore treated together. The optimal mea- 
surement angle depends on the measurement outcome, 
and this dependency is shown in Fig. [7] Because the op- 
timal measurement angle depends on the measurement 
outcome. Note that while all four measurement outcomes 
(there are two corresponding to s — 0) are equally proba- 
ble in the cluster state, the probability of the s = result 
vanishes for the GHZ state while the two remaining are 
equally probable. 

Alternatively, one can view the measurement as a weak 
measurement described by four POVMs described in 
terms of the three measurement angles, i9 p with p the 
measurement outcome, as shown in Fig. [7] The POVMs 
are thus 

E-i = cM_(i?_i) <g> M_(i?_ x ) 

E Q = cM_(tf ) <S>M + (tfo) 

££ = cM + (tf )®M_(tf ) (31) 

E 1 = cM + (i?i)<8>M+(i?i) 
Ex = 1 — E-i — E — E' — Ei, 

where M±{$) is a local projective measurement along 
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the direction and the constant approximately c < .5 is 
required for Ex to be positive in general. The element 
Ex , consisting of all local measurement outcomes that do 
not correspond to one of the others, is considered a 'fail- 
ure' outcome. If this measurement result is obtained, the 
result must be discarded, while for any other result, the 
measurement sequence disentangles the two outer qubits. 
The probability of failure is p(Ex) — (ip\Ex\ip)- For the 
specific model in question, c < .8 is sufficient to retain 
positivity for all J. Since a higher c means a lower prob- 
ability of failure, on can further minimize failure by se- 
lecting the highest c that allows for a positive Ex for a 
specific set of measurement angles. 

One may imagine that the POVM can be replaced by 
a projective two-qubit measurement, or equivalcntly that 
one can make adaptive measurements without postselec- 
tion. However, this is not possible for any < J < oo. 
Consider making two measurements in the four qubit 
chain with different angles, 9\ and 02- Mapping out 
the combinations that disentangle the chain for all four 
measurement outcomes, a pure projective measurement 
would be possible if all four lines crossed at some com- 
bination, while an adaptive or entangled measurement if 
one (say) 0\ for an outcome meant one could pick out a 
definite 02 depending on the outcome and be ensured dis- 
entanglement. However, our numerical investigation has 
determined that this situation does not occur for any fi- 
nite J. Hence, postselection is the only way to ensure 
disentanglement . 

Note that these results for the N = 4 chain cannot be 
immediately extrapolated to a general N qubit chain, but 
the same general results hold, and a disentangling mea- 
surement can still be constructed by equal measurements 
on the (N — 2)/2 middle qubit pairs with disentangling 
angles close to those described here. 
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APPENDIX A: JORDAN- WIGNER TRANSFORM 

The Jordan Wigner transform [34J is a mapping from 
the spins defined by the Pauli operators into spinless 
fermions, and is a standard technique in condensed mat- 
ter physics. We give a quick overview here, and show 
how to obtain the relevant correlation functions. For a 
N qubit chain, N spinless fermions are defined by anni- 
hilation operators 

and the adjoint creation operators. Futher one defines 
2N self adjoint Majorana fermions 



72^-1 



1 

iV2 



(a 



/< a l) 



1 

71 



where {7i,7j} = Sij. Switching coordinates x <-> z for 
the Hamiltonian ^ with B z = gives 



with the 2N x 2N matrix 



VII. CONCLUSIONS 
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For a ID spin chain for which the ground state is a 
cluster state, we have demonstrated the existence of a ro- 
bust "cluster phase" under local (B^ -field) and coupled 
(Ising) perturbations. All states in this phase exhibit di- 
verging localizable EL. However, as we have been unable 
to determine the bi-product unitary as a function of the 
measurement results, the usefulness of such states as a 
quantum computational wire may be limited. 

The existence of this cluster phase in one dimension 
may be extended to two dimensions even though our cur- 
rent numerical procedures are unlikely to be suitable for 
this. When extending to higher dimensions by similar 
schemes to MPS [32], the entanglement properties are 
likely to be even harder to distill. Alternative techniques 
beyond the MPS paradigm may be more suitable [33 . 



where 



B = 



-B a 
B x 



J = 



J 




1 = 



1 




and P = Yl eft = ± 1 is the parity of the model which is a 
conserved quantity, [H, P] — 0. If one considers a chain 
with open boundary conditions this simply amounts to 
setting P = 0. It is not immediately obvious which parity 
segment the ground state belongs to, but this is easy to 
verify by computing the two energies. For large N the 
difference between the two obtained states are in any 
respect very small. 

Define the correlation matrix Ty = ([7i,7j]), in 
which every second entry is zero since (72^721/) — 
(72^+172^+1) = 0. This enables us to compute the von 
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Neumann entropy of any subset of qubits effectively |35j , 
and expectation values also follow easily. 

The local z expectation value (which corresponds to x 
in the original coordinates) is simply (a*) = I^-i^ = 
Ti2 where the last step is valid under PBC, or P ^ 0. As 
an example of a more complicated expectation value, con- 
sider (16a), which in the current coordinates and trans- 



formed to Majorana fermions is 

° = < cr 2 cr 3 cr 5"- fJ W-l cr W> 

(N/2-2 
74fc+374(fe+l) ] 72JV)- 
fc=l 



Using Wick's theorem, and the fact that only odd/even 
combinations of indices couples, one finds that O = det G 
where G is a submatrix of T, 



G 



T3,8 
T7,8 
Til, 8 



^3,12 
IYl2 



T3.16 
IYl6 



T3,2Ar-4 
T7,2Af-4 



^3,2N \ 
^7,2N 



Vr 



2AT-5,8 T2Af-5,12 



r 



2Af-5,2AT 



Indeed, any expectation value is the determinant of a 
dense submatrix of T, the size of which depends on the 
expectation value in question. 
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